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' Non-gaussianity in the microwave background radiation is bound to play a key role in 

o ; . . 

^SJ , giving us clues about the physics of the very early universe. However, the associated calcu- 

Oh. lations, at second and even third order in perturbation theory, tend to be complicated to the 

' point of obscuring simple underlying physical processes. In this note, we present a simple 

analytic procedure for approximating the non-linearity parameters /atl and gisiL for cyclic 
models in which the cosmological perturbations are generated via the entropic mechanism. 
Our approach is quick, physically transparent and agrees well with the results of numerical 
calculations. 
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Observations of the cosmic microwave background radiation are quickly becoming detailed 



^ ' enough that within the next few years we can hope to obtain highly informative limits on the 

oo ; 

ly^ , bispectrum (and perhaps even the trispectrum) of primordial curvature perturbations [1]. In this 

in 

^SJ ' respect the detection/non-detection of non-gaussianity will provide a powerful tool in discriminat- 



ing between various theoretical models for the early universe. In simple inflationary models, the 
inflaton field is an almost free field, and correspondingly the curvature perturbations that these 
t> . models generate are governed by very nearly gaussian statistics [2]. More complicated inflationary 

^ ' models, such as multi-field models, can produce pretty much any value for the so-called "local" 

non-linearity parameters f^L ^-nd gj^i (corresponding to "squeezed" configurations in momentum 
space), which makes it difficult to predict a natural range [3]. So-called DBI models in which the 
inflaton possesses a non-canonical kinetic term lead to more distinct non-gaussian signals, involv- 
ing quite different ("equilateral") momentum configurations in their correlation functions [4], and 
multi-field DBI models can even produce significant contributions of both local and equilateral 
type simultaneously [5]. In ekpyrotic [6] and cyclic models [7], the cosmological perturbations are 
generated during a slowly contracting ekpyrotic phase with ultra-stiff equation of state w^k ^ 1 
(see [8] for a review). Such a phase can be modelled via scalar fields with steep negative poten- 
tials. The steepness of the potentials implies that these scalars are necessarily self-interacting, and 
this leads to natural values of the local non-linearity parameters that are in a range that will be 
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accessible to near-future observations [9-12]. Of course, this is why it is important to understand 
the physics that is responsible for the non-gaussian signals well, so that, in case of a detection, the 
consequences can be best appreciated. 

Here, we will focus exclusively on cyclic models in which the cosmological perturbations are 
generated by the entropic mechanism, as this is currently the best understood mechanism for 
producing a nearly scale-invariant spectrum of curvature perturbations during a contracting phase 
[13]. Recently, both the associated bispectrum and trispectrum have been calculated numerically 
[11, 12]. These calculations are rather involved, and do not provide many clues about the final 
outcome. Hence it is desirable to develop analytic methods, even though they might only be 
approximate, to understand the physics of these calculations more thoroughly. Some headway 
in this direction was made by the authors in a recent publication dealing with the bispectrum 
calculation [14]. In this note we present a new and much simplified approach that can in fact be 
applied to the calculation of the trispectrum as easily as to that of the bispectrum. 

The models under consideration can be described by gravity minimally coupled to two scalar 
fields with potentials. Nearly scale-invariant entropy perturbations are generated first, during a 
slowly contracting ekpyrotic phase that at the same time resolves the cosmological flatness puzzle. 
Subsequently, in the approach to the big crunch, the ekpyrotic potential becomes unimportant, 
and the universe enters a phase dominated by the kinetic energy of the scalar fields. During 
this phase, the entropy perturbations are converted into adiabatic curvature perturbations with 
the same spectrum, and these curvature perturbations form the seeds of the large-scale structure 
during the subsequent expanding phase. We will discuss the kinetic conversion phase in more detail 
below. 

We adopt the following parametrization of the potential during the ekpyrotic phase: 

Veu = -Foe^'^[l + es' + ^e^'^s^ + ^e's' + •••], (1) 

where we expect K3, K4 ^ 0(1) and where e ~ 0(10^) is related to the ekpyrotic equation of state 
Wek via e = 3(1 + Wf.^^jl. Wc use a to denote the adiabatic direction, i.e. the direction tangent 
to the scalar field space trajectory, and s to denote the "entropy" direction, i.e. the direction 
perpendicular to the background trajectory (note that the fields a and s arc thus defined such 
that the coordinate system they imply moves along with the background trajectory [15]); see Fig. 
1. The ekpyrotic potential is tachyonic in the entropy direction, and this instability causes the 
entropy perturbations to grow [13]. Moreover, this instability has the consequence that the global 
structure becomes a "phoenix" universe, in which the universe loses most of space to black holes 
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at the end of each cycle, while the regions that survive the big bang are aided by the dark energy 
to grow into vast new habitable regions - this was discussed in some detail in [16]. 

During the ekpyrotic phase, it is straightforward to solve for the entropy perturbation, with the 
result that [12] 

6s = SsL + S25sl + ssSsl, (2) 

with the linear, gaussian part Ssl being inversely proportional to time t (defined below in Eq. (9)); 
the coefficients S2 and S3 are given in terms of the parameters of the potential by 

S2 = — g— , (3) 

,3 = + ^ _ (4) 
^60 80 S'^ ^ ' 

The local non-linearity parameters /jvl and gjsfL can be defined via an expansion of the curvature 
perturbation ( in terms of its linear, gaussian part Cl, 

C = Cl + ^nlCI + ^9nlCI (5) 

Then, it was found numerically that, for conversions lasting on the order of one e-fold of contraction 
of the scale factor, the non-linearity parameters can be well fitted by the simple formulae [11, 12] 

3 

/atl ~ 12s2 + ^ = 2 '^^^ ^ 

The simplicity of the end result (given the complications of the third order perturbation equations 
involved) suggests that there ought to be a more straightforward way to obtain it. In fact, the 
physics of the kinetic phase (which follows the ekpyrotic phase, and during which the conversion 
takes place) is really quite simple, and moreover, except for the fact that its initial conditions involve 
the entropy perturbation Ss, the kinetic phase has no memory of the details of the ekpyrotic phase. 
In particular, only the total Ss in (2) matters, and the way we choose to decompose it into linear, 
second- and third-order parts is irrelevant at this point. This realization is the first ingredient of 
our calculation. 

The second is a compact and very useful expression for the evolution of the curvature pertur- 
bation ( on large scales and in comoving gauge [10, 17]: 

._JHSV_ 
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FIG. 1: After the ekpyrotic phase, the trajectory in scalar field space enters the kinetic phase and bends - 
this bending is described by the existence of an effective repulsive potential (the potentials are indicated by 
their contour lines). A trajectory adjacent to the background evolution can be characterized by the entropy 
perturbation Ss{tek-end) at the end of the ekpyrotic phase, leading to a corresponding off-set 6s{thend), or 
equivalently SV (tbend) , at the time of bending. 

where a dot denotes a derivative w.r.t. time t, H = a/ a is the Hubble parameter and a the scale 
factor, 5V = V{t,x'^) — V{t) and a bar denotes a background quantity. This equation is exact in 
the limit where spatial gradients can be neglected, and can thus be expanded up to the desired 
order in perturbation theory if required. First, let us present its derivation [17]: considering only 
very large scales, we can write the metric as 



where all the inhomogeneities are in C,. This defines Q to all orders in the long- wavelength limit. 
Then the equation of continuity reads 



where p is the (scalar) matter energy density, and P its pressure. But we're interested in the 
curvature perturbation on surfaces of uniform energy density, so p = p (and hence also H = H). 
And since p satisfies p + 3H{p + P) = 0, we immediately obtain 



(9) 



p + 3{H + C){p + P) = 0, 



(10) 



C = -H 



5P 



(11) 



P + P + 5P 
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Now, since we choose to consider hypersurfaces on which 5p = 0, we obtain the relations 5{a'^) = 
—2dV and thus SP = —2SV. Plugging these relations into (11) then yields our desired result, Eq. 
(8). 

Incidentally, by expanding Eq. (8), it is also possible to show that it is equivalent to the third- 
order equation derived in [12] using the covariant formalism. The crucial thing is to keep in mind 
the definitions of the adiabatic perturbation 5a and the entropic one 6s at higher orders, provided 
in [15] and [12]. In particular, we have (from hereon we drop the bar on background quantities - 
this should not lead to any confusion) 

SV = Vs6s + VJa 

+Vs[Ss^^^ + ^SaSs + ^i6af] + V^iSa^''^ - ^SsSs] 
' a 2a ' 2a 

+ ■■■ (12) 
Up to third order and in comoving gauge {Sa = (5(7(2) ^ ^^(3) ^ 0) we then get 
SV = Vs6s 

+y^5,{3) _ v^[^{SsSs^-')y+ ^{SsfSs] + VssSsSs^''^ - ^{SsfSs + ^F,,,(<5s)^(13) 

and by expanding (8) we obtain 
2 ff 

C = -tVsSs 

+ ^-^[V'js6s('^ - ^V,VASs)'6s + Iv^VMSs)'] + ^-^V^siSsf, (M) 

which agrees precisely with the equation derived in [12]. Having shown this equivalence, we will 
now stick with the simple and compact form (8). 

We are assuming that, during the kinetic phase, the trajectory in scalar field space contains a 
bend, and this bend is what causes the entropy perturbations to source the curvature perturbations. 
In cyclic models embedded into heterotic M-theory, such a bend occurs naturally because the scalar 
field space contains a boundary which effectively acts as a repulsive potential (we refer the reader 
to Ref. [18] for details). However, we are simply citing this example as a concrete realization. 
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Our calculation applies to all cases where there is a bend in the trajectory, although for extreme 
(and unnatural) cases where the bending angle is close to 0° or 180° some of our approximations 
below might break down. We are assuming that this bend can be described as being caused by a 
monotonic repulsive potential, as depicted in Fig. 1. 

Now, the third and last ingredient of our calculation is the simple relationship between 6V and 
5s during the conversion process. During the ekpyrotic phase, the curvature perturbation picks 
up a blue spectrum [19] and is hence completely negligible on large scales. To be precise, since 
SV 7^ during ekpyrosis, there is already some conversion of entropy into curvature perturbations 
occurring at this stage. However, this contribution is entirely negligible compared to the subsequent 
conversion (see [12, 14] and note that since = during ekpyrosis, SV starts out at subleading 
order), and hence we can take C(^efe-end) ~ where tek-end denotes the time at the end of the 
ekpyrotic phase, or equivalently, at the start of the kinetic phase. Moreover, as we will see below, 
at the end of the conversion process (" is still significantly smaller than 6s, and hence, during 
the conversion process, we can take the potential to depend only on 6s. And since the repulsive 
potential is monotonic, and we are interested in small departures Ss <^ 1 from the background 
trajectory, it is intuitively clear that SV is directly proportional to Ss during the bending. A 
numerical calculation readily confirms this simple relationship. 

During the conversion, the effect of the repulsive potential is to cause the entropy perturbation 
to behave approximately sinusoidally, independently of the precise functional form of the potential 
(this was shown analytically in [14]): 



where tc denotes the time at which the conversion starts. Moreover, the precise value of Ss{tc) is 
unimportant for the present calculation. The frequency lo is of 0(1/ At), where At is the duration 
of the conversion; the more careful analysis presented in [14] leads to u> ^ 2.5/ At. Another useful 
quantity is the rate of change of the angle of the trajectory in scalar field space [20] 

^ = _^~ 1. (16) 
& At ^ ' 

Also, the scale factor and the scalar field velocity along the background trajectory are rather 
unaffected by the presence of the repulsive potential, so that they simply assume the values they 
would in the absence of any potential 



Ss cos[Lo{t — tc)]Ss{tc), 



(15) 




-V2 



(17) 



a = 



Vst ■ 
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As we will confirm below, during the conversion process SV <^ so that Eq. (8) simplifies 
further to 

2 ff 

C ~ JT^V. (18) 
Then, at linear order, we immediately obtain 

Cl = / -^e6sL (19) 

J bend ^ 

sin(a;Ai)(5s(tc) (20) 

3 u) 

~ \Htc). (21) 
5 

But, as argued above, 5s as a whole must behave approximately in this way during the conversion 
phase, and subsequently analogous relationships hold at higher orders too: 

C(^) « \s25sl C^^) ^ \s,5sl. (22) 



These expressions immediately allow us to calculate the non-linearity parameters 

5 C(2) , 

Inl = ;^"72" ~ §552 ~ 8s2 (23) 
25 

g^^ = —^ ^ f 52^3 « 70S3. (24) 



Thus, without much work at all, and to better accuracy than a factor of 2, we recover the numer- 
ically (and laboriously) obtained fitting formulae in Eqs. (6)-(7) above. 

Before discussing this result, let us briefly pause to verify the approximation made in obtaining 
Eq. (18): during the kinetic phase, we can rewrite (8) as 

C-T^- (^^) 

The approximation made above consists in writing ( ^ tSV and this leads to C ~ ^'tlend^^ i^bend) ■ 
But we know that by the end of the conversion process ^ and hence we find that 

^tlndSV ^ 6s « 1, (26) 

which shows that the approximation is self-consistent and confirms the validity of (18). 

So what does our result tell us? The main point is that due to the simplicity of the kinetic phase, 
the non-linearity that was present in the entropy perturbation gets transferred straightforwardly 
to the non-linearity in the curvature perturbation. Our calculation therefore explains why there 
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are no significant additional constant terms in (6) or constants and Ks-dependent terms in (7); a 
priori, there was no reason for such terms to be absent. 

Moreover, the overall magnitude of f^L and qnl '^^ set solely by the efficiency of the conversion 
process, as expressed by the relationship between 5sl and Cl in Eq. (21). The fact that no 
additional parameter enters into Eq. (21) has important consequences for observations, as it 
determines the scaling of the non-linearity parameters with the equation of state parameter e, as 
expressed in Eqs. (6) and (7). A natural value for e would be about 50, so that we can expect 
/iVL to be of order a few tens, with the sign typically determined by the sign of K3, and ^atl to 
be of order a few thousand and typically negative in sign. These values represent the natural 
values predicted by models making use of the entropic mechanism. They comfortably fit current 
observational bounds [21, 22] while being detectable by near-future observations. It is useful to 
contrast these values with those predicted by "new ekpyrotic" models [23, 24] where the entropy 
perturbations are converted into curvature perturbations directly during the ekpyrotic phase. For 
this variant conversion process, the dependence on the equation of state e is more pronounced, 
with /jvL oc e and g^L oc [9, 10, 12]. In addition, f^L is predicted to take a negative value and 
qnl ^ positive one; the magnitude and sign do not fit well with current observations. 

Finally, we note that, in spirit, our approach is somewhat reminiscent of the 5N formalism 
[25, 26]. However, a full 6N calculation spanning both the generation and conversion of the 
cosmological perturbations is made difficult here because of the transition between the ekpyrotic 
and kinetic phases; note, in particular, that in going from the ekpyrotic to the kinetic phase, the 
equation of state drops drastically from Wf^k 3> 1 to Wkin ~ 1 (by contrast, the 5N formalism is 
well adapted to new ekpyrotic models where the generation and conversion both take place during 
the ekpyrotic phase [9]). In fact, it is precisely the disconnectedness between the two phases that 
allows our method to work so well, and we would expect it to be applicable more generally to 
cases where the physical properties of the phases of generation and conversion differ substantially. 
The separation between the two phases allows for a two-stage approach in which we first solve for 
the entropy perturbation during the ekpyrotic phase, and then use this as input for calculating 
Q or, cquivalcntly 5N, by perturbing around the background trajectory during the kinetic phase. 
The simple, yet non-perturbative, Eq. (8) then reveals its full effectiveness by yielding the result 
in just a few lines of derivation. In this way, we have found a quick and rather accurate way of 
understanding non-gaussianity in two-field cyclic models of the universe. 
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